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BUMPER-PROTECTED LAMINATED SPACECRAFT MAINWALLS 


By Ray Kin slow 

Department of Engineering Science 
Tennessee Technological University 

1. INTRODUCTION 


It was first suggested by Whipple (Reference 1) that a thin shield spaced 
some distance from the mainwall of a spacecraft would completely fragment a 
meteoroid and disperse the resulting debris over a large area of the main struc- 
ture inflicting only slight damage. Many subsequent investigators have con- 
firmed the fact that such a shield or bumper can significantly reduce meteoroid 
impact damage to the structure. 

It has also been demonstrated that laminated mainwalls can be utilized to 
reduce the damage resulting from hypervelocity impact (Reference 2). 

This study is an investigation of the effects of employing the use of both 
a bumper and a laminated spacecraft mainwall to reduce the meteoroid hazard to 
space vehicles. A quasi- theoretical approach has been taken. A mathematical 
model has been formulated that gives results that are in good agreement with 
experimental data. 

2. DEBRIS CLOUD PRODUCED BY PROJECTILE-BUMPER INTERACTION 


As previously stated, the bumper serves to fracture the hypervelocity 
projectile and to diverge its fragments. Various assumptions and computations 
as to the debris cloud characteristics have been made by many different 
investigators . 

It is assumed that the origin of the debris cloud is one projectile 
diameter in front of the bumper. This is the conclusion reached by several 
investigators. Turpin (Reference 3, p. 118) says, "This location appears to 
vary slightly with material, but generally corresponds to the rear surface 
of the projectile at the first instant of impact." According to Rosenblatt 
(Reference 4, p. 69), "The velocity distribution in the debris cloud indicates 
that the mass elements making up that cloud have the appearance of emanating 
from just ahead of the target surface. The average intercept of the velocity 
vectors with the axis is about one projectile diameter ahead of the target for 
all impact cases considered." This assumed point .of origin of the debris is 
indicated by point O' in Figure 1. Its distance from the main structure is Y, 
the sum of bumper spacing, bumper thickness, and projectile diameter. 

In order to treat the problem theoretically, Maiden, et al (Reference 5, 
p. 21) assumed that the fragments loaded the mainwall uniformly over a circular 
area of diameter equal to one-half the bumper spacing. This is indicated by 
the rectangle in Figure 2 which shows the normalized axial component of momen- 
tum of the debris per unit area. In a later study of the ballistic limit of 



double-walled meteoroid systems. Madden (Reference 6, p. 9) assumed that the 
mass distribution of the spray particles could be approximated by a normal 
Gaussian distribution. Assuming that the momentum applied to the mainwall has 
the same distribution as that of the mass discussed by Madden, it can be de- 
scribed by the relation, mom = This is also shown in Figure 2. 

In an analytical study of debris clouds due to impacts on thin plates , 
Rosenblatt, et al utilized the STEEP code to obtain numerical solutions (Ref- 
erence 4). The axial momentum per unit area produced by the impact of a 1/8- 
inch aluminum sphere on a 1/32-inch aluminum plate at a velocity of 24,600 ft/ 
sec is shown in Figure 2. This is in close agreement with the Gaussian distri- 
bution, especially for larger angles of dispersion. 

It was pointed out in Reference 4, page 35, that three regions of the 
momentum (as well as mass and energy) distribution curves can generally be 
recognized. The first region is due to the main body of the debris cloud and 
the most severe damage to the mainwall will be caused by the impact of these 
fragments . This region is confined to an included half-angle of less than 
20 degrees . 

The momentum distribution per unit area used in this study is assumed to 

be proportional to (1 - A tan 2 0), where 0 is the included half-angle and A is 
a debris distribution factor. This is shown in Figure 1 with A having a value 
of 6. This is in good agreement with the values computed for the Al-Al impact 
for angles up to about 17 degrees which includes the area where damage is most 
likely to occur. 

This more simple expression was used because: (1) For the material com- 

bination being considered for the laminated mainwall, there is a total reflec- 
tion of the longitudinal waves at an angle of about 25 degrees; (2) It has been 
pointed out that damage will probably be limited to an included half-angle of 
less than 20 degrees; and (3) Mathematical and computer difficulties were en- 
countered in assuming a Gaussian distribution. 

It is assumed that the leading edge of the debris cloud travels at the 
projectile impact velocity. According to Reference 4, page 58, the ratio of 
cloud velocity to projectile velocity varies from 0.92 to 1.0 for Al-Al impact 
at velocities of 0.75 and 1.5 cm/ysec, respectively. This ratio is slightly 
greater than unity for copper and cadmium impacts. Swift and Prater (Refer- 
ence 7, p. 11) give a cloud- to-pro j ectile velocity ratio of about 0.95 on the 
cone axis to 0.85 at an angle- of 20 degrees for Cd-Cd impact at 0.67 cm/ysec. 
This assumption appears to be reasonable in view of both experimental and 
theoretical results . 

Many excellent photographs have been made of debris clouds. Figure 3, 
taken from page 3 of Reference 8 is one such photograph-. As the envelope of 
the debris resembles a Lemniscate of Bernoulli, it was suggested that this 
curve might be used as a mathematical model of the cloud. True Lemni.s cates 
are shown in Figure 4. Although these resemble the cloud photographs, the 
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angle of 90 degrees formed by the tangents at point 0 is too large. Several 

modifications are shown in Figure 5. The curve, r i 2 = cos38, most accurately 
matches the cloud photograph. Figure 6 shows the cloud profiles at various 
times based upon this model. 


If the projectile and debris velocity is denoted by v and time is measured 
-fr om -the instant of the first mainwall impact, the time at which any portion of 
the cloud at a cone angle of 0 reaches the wall is given by the relation 


t 


Y 

v 


r i f 

LcosO ( cosBQ ) n 


3 


( 1 ) 


where B and n are spray dispersion factors, having values of 3 and 1/2, re- 
spectively, in this example. 


3. MAINWALL LOADING 


The loading of the spacecraft's hull by the fragments in the debris cloud 
is a function of many factors including the mass and velocity of the projec- 
tile, the bumper spacing, the debris distribution, the cloud dispersion, and 
the momentum multiplication. 

The characteristics of the debris cloud have been discussed. The momen- 
tum multiplication factor results from two multiplication effects . The first 
results from the projectile impact on the bumper and the second results from 
the reflection of gas and debris from the hull. Cunningham (Reference 8, p.47). 
Maiden (Reference 5, p. 21), and many other investigators agree that for most 
metals the momentum transferred to the mainwall is equal to approximately twice 
the momentum of the original projectile. 

Because the debris cloud is practically hollow and its wall thickness is 
small (Reference 4, p. 57), and also because the fragments are very small, it 
is obvious that the loading time at any point is extremely short (Reference 5, 

p. 21). 

It is assumed that the mainwall pressure per unit area can be specified 
by the relation 


i \* -a t -a t 
p -v e i L - e 2 

where and a are decay constants and t is time; For any desired rise and 


( 2 ) 


decay times , the values of the decay constants can b.e determined .from the chart 
of Figure 7. 

The normal pressure (see Figure 1) on the mainwall surface is 

p = K( 1-A tan 2 0)(e _cx l t - e _ 01 2't) ( 3 ) 
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( 4 ) 


Py = K cos0 (l-A tan 2 0)(e -a l t - e a 2 t ) 
where K is a constant depending upon several parameters. 

Equating the product of the projectile momentum and the momentum multi- 
plication factor to the total normal component of the momentum applied to the 
mainwall gives the relation 

M(mv) = /“ /tan-V 0 ^ da dt (5) 

t=o 0=0 y 


where 


da = 2nv dr = 2 ttY 2 sec 2 0 tan0 d0 
x x 


therefore 

M(mv) = 2uKY 2 //(l-A tan 2 0)(e - e -a 2 t ) sec0 tan0 d0 dt 


= 2irKY 2 


r a 2 -ot ii 

■ 




2(A+1) 1 - 5 3+2A' 


3A0.5 


If the debris distribution factor is 6 (max 0 = 22.2°) 

M(mv) = 0.307 KY 2 


V“i 


a l a 2 


( 6 ) 

(7) 

( 8 ) 

(9) 


4. PROPAGATION OF SPHERICAL STRESS WAVES 


The debris particles impacting the mainwall surface generate spherical 
dilatational stress waves in the structure. This wave propagation can be 
specified by the equation 


a 2 <j> 2 /ii± 2 a$\ 

3t 2 " Cl Ur 2 r 3r/ 


( 10 ) 


where is a scalar displacement potential, c^ is the wave velocity, r is the 

radius vector from the point of impact, and t is time. Particle displacement 
(u) is specified by the relation 


u 


34 >_ 

9r 


( 11 ) 


The radial and tangential stresses are 
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(12) 


a 

r 


(A+2y ) +2A 



and 


a 


0 


A 





(13) 


where A and y are the Lame constants and are related to Young's modulus (E) 
and Poisson's ratio (v) as follows: 


(ltv)(l-2v) (14) 

u = ( 1 5 ) 

If the time- varying pressure (p) given by Equation (3) is applied to the 
surface of an assumed hollow hemispherical cavity of radius R q in the target 

with its center at the point of impact, elastic stress waves will be gener- 
ated in the structure. This mathematical model is that described by the 
author in Reference 9. 


A solution of the wave equations, based on Blake's work, (References 10 
and 11) is ’ 


* = l . ta £ j)ro I t+ e - aoT h + Krfy\ 2 ' 

i=lP r ( w o + ( a o _a i ) 2 L V “o / . 


where 


0.5 


cos 


0 ) o t- tan" a ^' 
0)0 / 


(16) 


J J 


t = t - \ 

“l I 


_ C 1 / 1- 2v\ 

“o " ^ 


w o = 


_ C 1 


(l-2v) 


0.5-, 


1-v 


(17) 

(18) 
(19) 


5. REFLECTION AND REFRACTION AT AN INTERFACE BETWEEN TWO MEDIA 


When an elastic wave of dilatation reaches a boundary between two media 
as is the case of the laminated mainwall considered in this study, four waves 
are generated. Two of these, a dilatation or longitudinal wave (RL) and a 
distortional or shear wave (RS), are reflected back into the first medium, 
and two such waves (TL) and (TS) are transmitted into the second medium (Ref- 
erence 12). 
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At the interface there are four boundary conditions. These are that on 
both sides of the boundary the following quantities must be equal: 

(a) The normal displacements; 

(b) The tangential displacements; 

(c) The normal stress; 

(d) The tangential stress. 


The transmitted and reflected waves and their angles with a normal to the 
interface are shown in Figure 8. RL, RS, TL , TS are the stress amplitudes of 
the waves . 

The relations between the wave velocities and the angles of incidence and 
reflection are: 


sin 0 sin 3> sin B sin 0' sin 3' 


( 20 ) 


"1 "1 1 2 2 

where c^ and c^ are the velocities of the longitudinal and shear waves in the 


first medium, and c 2 and c£ are the corresponding velocities in the second 
medium. These velocities are given by the relation 


C 1 = 


= 


Co = 


Ei(l-vi) 


jp 1 (l+v 1 )(l-2v 1 )J 
-,0.5 


0.5 


Co = 


2p 1 (l+v 1 ) 


~ E 9 (1-v 9 ) 

l P 2 (1+v 2 )(1-2v 2 ) 

r Eq 1 °’ 5 
2p 0 (l+v o ) 


0.5 


L M 2 


2 J 


( 21 ) 

( 22 ) 

(23) 

(24) 


The four boundary conditions lead to the following four equations which, 
together with the relations of Equations (20), may be solved for the stress 
wave amplitudes and their directions. 


t 


cos 8 
p l c l 


RL 

P 


sin B 

_ RS n 

cos 0' 

'tl' 

p i c i 

_P J 

* P 2 c 2 

P . 


sin B' TS~ j 

+ P 2 C 2 L.P J 


cos 0 

+ Pl c l 


(25) 


sin 0 

"rl" 

j cos B 

i 

CO 

04 

sin0 ' 

’tl' 

cos B ' 

fTS 

+ p l c l 

Lp . 

! + Pl c^ 

LP 

P2 C 2 

lP J 

P2 C 2 ! 

1 

LP _ 


- cos2B 



+ sin2B 



+ cos2B' ! 


fTL 


+ sin2B’ 


fTS] 

Lp j 


sin 0 
p l c l 

+ cos2B 


(26) 

(27) 
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,! \2 

+ (— i- \ sin20 
lc l 


RL 


+ cos2gW + MV sin2e* 


IP J 


'- c 2 


i_p j 


- cos2B'i 


TS 


= 


J 


-li 


sin20 (28) 


The "two materials being considered for laminated spacecraft mainwalls are 
polycarbonate and aluminum. Nominal values of the mechanical properties of 
these materials are: 


Aluminum 

E-= 6.9 x 10 11 dynes/cm 2 

p = 2.68 

v = 0 . 33 


Polycarbonate 
E = 0.24 x 10 11 dynes /cm 2 

p = 1.21 

v = 0.40 


The outer layer of the material is polycarbonate and the material on the 
inside is aluminum. The incident stress wave at the interface is assumed to 
have an amplitude of unity. The values of the reflected and transmitted wave 
amplitudes and the angles of reflection and refraction are given in Figure 9. 
When the velocity of propagation of the transmitted wave is greater than that 
of the incident wave, as in this example, there is a critical angle of inci- 
dence which makes the angle of refraction 90 degrees . This angle is approxi- 
mately 22 degrees for the polycarbonate- aluminum laminate being considered. ■ 
For angles of incidence greater than this, the relations that have been derived 
do not apply. For angles of incidence greater than the critical angle, the 
problem must be treated in terms of complex quantities (Reference 13, p. 38). 
This has not been done in this study. 

If the materials are reversed so that the outer layer is aluminum and the 
inner material is polycarbonate, the transmitted stresses are greatly reduced. 
This has been demonstrated experimentally by Gehring (Reference 14, p. 474). 
Figure 10 gives the results for this condition. In this case there is no 
critical angle and therefore no limitation on the value of the incident angle. 
Figure 11 shows the transmitted and reflected stress amplitudes as vectors, 
giving a clearer picture of both and angles and amplitudes. 

Although these results give a graphical comparison of the wave amplitudes 
they do not present the true picture when the stress waves result from the im- 
pact of fragments in the spray cloud. Assuming that the momentum distribution 

per unit area is proportional to (l-6tan 2 0) as shown in Figure 2, the values 
of the incident, transmitted, and reflected waveq will be as shown in Figure 
12. These values are shown as vectors in Figure 13. 

The location of the wave fronts in a polycarbonate- aluminum rnainwall at 
times from 0 to 2.0 microseconds are shown in Figure 14 for the conditions 
indicated. This figure is for illustrative purposes only and does not repre- 
sent actual dimensions. 
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6. METHOD OF STRESS ANALYSIS 


In Figures 15» and 16, Y is the distance from the- apparent point of origin 
of the debris fragments from the mainwall which is assumed to be equal to the 
sum of the bumper spacing, bumper thickness , and projectile diameter. The 
thickness of the first layer is indicated by T^ and the thickness of the second 

layer by The location of any point (P) in the material is given by the 

coordinates x and y. The y-distance is measured from the material interface. 

A positive value of y indicates that the point is in the second layer of the 
material and a negative value indicates that it is in the first layer. For 
the present, only points in the second layer are being considered. 

As fragments 'impact the mainwall it is assumed that spherical stress waves 
are generated but only those portions in the direction of the impacting parti- 
cles which will ultimately reach point P are being considered. There are, of 
course, an extremely large number of stress waves in the structure, one being 
generated by each debris particle. These will, according to the Huygens prin- 
ciple, result in essentially a single wave front as shown in Figure 14. It 
is assumed that the material at point P will be affected by only those por- 
tions of the waves that pass through the point. A brief outline of the methods 
used to compute the principal stresses within the second laminate will first be 
described, followed by the relations and equations used in the computations. 

The fragments moving in the direction 0 will generate a wave having an ampli- 
tude of p 1 when it reaches the interface. Two waves, a longitudinal and a 

shear, will be reflected at the interface and two such waves will be transmit- 
ted. The longitudinal wave (TL)^ will reach the point as shown in Figure 16. 

As was seen in Figure 13, the amplitude of the shear waves will be small com- 
pared with the longitudinal waves and their effects will be neglected. The 
fragments moving in the direction 0^ will likewise generate both shear and 

longitudinal waves at the interface. Considering only the longitudinal wave, 
it will be reflected from the inner surface to point P. It has been shown 
(Reference 15) that the shear wave reflected from a free surface will contrib- 
ute very little to the total stress. Each of the waves (TL)^ and (TRL)^ 

creates radial stresses (a ) and tangential stresses (a ) at the point. The 

r 0 

x and y components of stress due to each wave are computed and their combined 
effects are determined by superposition. From this, the maximum and minimum 
principal stresses and their directions are computed. This will give the 
stress history at that point. By considering a large number of points within 
the material, the stress field can be plotted for any time, as well as the 
maximum tensile and compressive stresses that are generated within the material. 

In this analysis, the forcing function is applied to the surface of a 
hemispherical cavity , the radius of which is equal to or greater than the dis- 
tance at which the stress wave may be considered to be elastic. The method 
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of determining the radius of this cavity (R^) will be discussed later. 

To simplify computations, it is assumed that the stress waves originate 
and travel through a homogeneous target of the second material. The imaginary 
cavities in which the waves originate are located so that the waves reach 
point P at the same time and from the same direction as if they moved along 
the paths shown in Figure 16. 


The center of the first cavity is located at 0, at a distance R from 
point P. 1 1 


R i = °i p = ^7 


COS0- 


. _ y 

r COS0j^ 


(29) 


at an angle of 0^ (Figure 17). Its radius is (c^/c-^. This wave will have 
the correct amplitude if the factor (TL)^/p^ is applied. Likewise, the second 
wave may be assumed to originate in a cavity of radius (c 2 R 0 /c 1 ) with its cen- 
ter (Og) at a distance from point P. 

c o r T n ~i r2T 0 -vi 

(30) 


r 2 = °2 P = c 


C 2 

Tl 


2T 2 -y 

= C 1 , 

cos 0 2 

t 

cos ©2 


at an angle of 9), • The correction factor for the wave amplitude is (TRL) 2 /p 2 , 

The following values must also be computed in order to calculate the 
principal stresses at point P. 


Values of 9^ and 0^ can be found by solving the two equations 

X = (YtT^) tan 0^ + y tan 0^ 


and 


'i 


and 


c l/c 2 = sin 0^/sin 0J 
Values of and 0), can be computed by solving the two relations 

X = (Y+T-^) tan 0 2 + (2T 2 ~y) tan 0), 

c l^ c 2 = s ^ n ®2^ s ^ n S 2 

Values of 8^, 8^, 6 2 » $ 2 9 an< ^ ^2 a:re determined from the relations 

sin 6 1 /c 1 = sin = sin B|/c^ 


(31) 


(32) 


(33) 


(34) 


(35) 
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(36) 


sin e 2 /fc l = sin = Sin ^V°2 


and 


sin 0^/c 2 = sin B^/c^ 

The times t and at which the forcing function is applied to the 
cavities 0 and 0 2 are 


(37) 


Y+R„ 


t l = 


cos0^(cos3e^ ) 


0.5 


- 1 


and 


Y+R. 


t 2 = 


,0.5 


- 1 


cos02 (cos302 y 

The times of arrival of the two wave fronts at point P are 


t, + 


T l- R o Jt 2-y 


u 1 r c. i i T 

- Lr - 1 Cj_cos0j C2Cos0-[ 

t = t + T i" R ° 2T 2- y , 

2f 2 CtCOsOq + CoCosBo 


(38) 


(39) 


(40) 

(41) 


Stress wave amplitudes (TL)^/p 1 and (TL) 2 /p 2 can be found by solving 

Equations 25, 26, 27, and 28 when the corresponding subscripts are used. 

The ratio of the incident to reflected waves (TRL) 2 /(TL ^ at the inner 

surface is found by applying the conditions of zero normal and shear stresses 
at a free surface. This gives the relation 

(TRL) 2 - sin20 2 sin202 + (c 2 / c 2 ) 2 cos 2 2g 2 ^ 

(TL) 2 + sin20 2 sin2B 2 + (c 2 / ’c^) 2 cos 2 2f3 2 

The ratio (TRL^/pg may now be found 


(trl) 2 (trl) 2 t (tl) 2 
p 2 (tl) 2 p 2 

The pressure pulse to be applied to cavity 0^ is that of Equation 
plied by the factor (TL)^/p^ giving 


(43) 

(3) multi- 
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K(l-A tan 2 0 1 ) e~ a l (t t l ) - e a 2 (t t l ) 


and the pulse to be applied to 0 ^ is 




-==== — K(d;-A— tan— 0 2 7 )— |e— 

p 2 - L 

__ __ 


The radial and tangential components of stress at point P produced by 
these waves are cr , a*. , o^ , and o + . In order to determine the resulting 
r i t i r 2 t 2 

stress it is necessary to resolve these stresses into their x and y components 
as follows 

°r 1 +° 0 1 a r 1 -a e 1 

a x . = o o cos 20 I 

X 1 2 ^ 1 


c r 1 +o 0 1 °r 1 _a 0 1 


cos 20 ! 


a _o ft 

r l 0 1 

_ _i _ — — s i n 2 0 1 
x-Yl 2 1 


= a r 2 to e 2 “r 2 -“e 2 


Ov, -0( 


cos 20 ^ 


a r 2 " a 02 a r>2 _a 02 

Oy 2 = 2 + j cos 20 £ 

a r 2 - a 0 2 

T = — t £ s in2 0' 

T xy 2 2 2 


The resulting x and y components of stress at the point are then 

o=o to 

X X-^ X 2 


0=0 to 

y yi y 2 


T =- T + T 

xy l xy 1 L xy 2 


From which the principal (p and q) stresses are found 


a x +a y 

P =- 2 “ + 


-i 0.5 


+ (x ) 2 

xy 
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1. Specify values of the constants E^» E 2 , pj_, p 2 , v-p V 2 , m, v, M, Y, T^, T 2 , 
R Q , A, B, n, a^j <* 2 ’ At, x, and y. 

2. Compute K - Equation (8). 

3. Stress wave velocities c, , c’ , c , and c' - Equations (21), (22), (23), 
(24). 1122 

4. Compute angles 0^ and 0| - Equations (31), (32). 

5. Compute angles 0 2 and 0 2 - Equations (33), (34). 

6. Compute angles (3^ and - Equation (35). 

7. Compute angles 8 2 and g 2 - Equation (36). 

8. Compute angle B 2 - Equation (37). 

9. Compute the time t^ - Equation (38). 

10. Compute the time X - Equation (39). 

11. Compute the time t^ - Equation (40). 

12. Compute the time - Equation (41). 

13. Compute the cavity radius = R o2 = c 2 R Q /Cj_ (page 9). 

14. Compute the distance R^ - Equation (29). 

15. Compute the distance R^ - Equation (30). 

16. Compute the ratios (TL) /p and (TL) /p - Equations (25), (26), (27), 

(28). 11 22 

17. Compute the ratio (TRL^/dL^ - Equation (42). 

18. Compute the ratio (TRL^/pg - Equation (43). 


12 



19. Compute the forcing function - Equation (44). 

20. Compute the stresses a and o A at time increments of At starting at 

r l °1 

time t-j^ - Equations (12), (13). 

21. Compute the x, y components ° x j 0 y » and T xv - Equations (46), (47). 


(48). 


x Yl 


22. Compute the forcing function F2 - Equation (45). 

23. Compute the stresses a and at time increments of At starting at 

2 2 

time - Equations (12, (13). 

NOTE: The stresses of step 20 must be computed for the same times of 

step 23 after time 

24. Compute the x, y components o x , o ' , and t - Equations (49), (50), 

(51).' ^ ^ ^ 

25. Compute the values of a , a , and t - Equations (52), (53), (54). 

X V 

26. Compute values of the principal stresses, p and q - Equations (55), (56). 

27. Compute the direction of the maximum principal stress, y - Equation (57). 


8. EXPERIMENTAL DATA 


As the method employed in this study is semi-empirical or quasi- 
theoretical in nature, some experimental data must be used if actual values 
of stress rather than relative values are to be computed. This experimental 
information was supplied by the NASA Langley Research Center (Reference 16). 
The information was for six targets, four of which were solid 2024-T3 alumin- 
um targets and two of which were polycarbonate- aluminum laminated targets. 

The data supplied are given in Table I . 

A search of the literature failed to reveal any significant information 
about the dynamic strength of this material. The peak tensile stress for 
incipient spall failure in 6061-T6 aluminum was given as 15 kilobars for a 
pulse duration of 0.1 microsecond in Reference 17, p. 398. The static ten- 
sile strength of 2024-T3 is about 50 percent greater than that of 6061- T6; 
however, this is perhaps somewhat greater than the dynamic strength ratio of 
these materials. It is believed that a dynamic tensile strength of 20 kilo- 
bars or about 300,000 psi is a realistic value and will be used in this study. 

Although a stress criteria for fracture is assumed in this analysis, it 
is realized that fracture is probably time- dependent and that the strain-rate 
would perhaps be a more significant factor to consider. 
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9. TENSILE STRESSES IN SOLID ALUMINUM TARGETS 


Because of the small size of the individual particles in the debris cloud, 
it is believed that the shock waves generated will decay into elastic waves 
within a very short distance. Because of this, it is concluded that the value 
of R q should be small, perhaps a fraction of a millimeter. In this study the 

value of R will be varied from 0.25 to 0.50 mm in order to find the value that 
o 

will give results that are in agreement with the experiments . 

Another factor that will affect the values of the computed stress is the 
shape of the forcing function applied to the cavity radius . The forward momen- 
tum of this forcing function will be equal to that of the impacting debris but 
its amplitude will depend upon the rise and decay times of the pulse. Rise 
times of 0.05 and 0.10 microsecond will be investigated and the ratio of the 
decay to the rise time will be assumed to have a value of five. 

The method given for the analysis of stresses in laminated targets may be 
used for solid targets by assigning a value of zero to T . 

The target dimensions are those given in Table I for shots 1, 2, 3, and 
4. The value of Y is 1.00 cm, the sum of the projectile diameter (0.153 cm), 
the bumper thickness (0.043 cm), and the bumper spacing (0.800 cm). As the 
projectile velocities are approximately the same, an average of 7.85 km/sec 
is used. 

There are three possible causes of tensile stress, any of which may be 
the maximum. These are illustrated in Figure 18. There will be tangential 
tensile stresses near the front surface of the target. The radial fractures 
seen extending out from the crater in transparent plastic targets are produced 
by these stresses. There will be radial tensile stresses produced by the tail 
of the forward-moving wave. These are usually of relatively small amplitude. 
There will be the radial tensile stresses caused by the reflected compression 
wave. These stresses are the cause of internal cracks approximately parallel 
to the rear surface or the spallation of the target's rear surface. The max- 
imum of these tensile stresses is always considered at any point. 

Figure 19 shows the maximum tensile stress along the y-axis of the four 
solid aluminum targets using a forcing function of 0.05-microsecond rise time. 
The effect of the value selected for R^ is apparent. If R q is assumed to be 

0.050 cm, the material must have a dynamic tensile strength of approximately 

114 kilobars as shown because target No. 3 failed by spallation and target 

No. 4 did not fail. Likewise, for a value of R of 0.0375, the tensile 

o 

strength must be about 80 kilobars; and for R q = 0.025, the strength must be 
about 50 kilobars. 
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Changing the decay constants so that the forcing function has a rise time 
of 0.10 microsecond, a value of 0.050 cm for R Q would require a dynamic tensile 

strength of about 54 kilobars. Decreasing the cavity radius to 0.0375 cm de- 
creases the required strength to about 37 kilobars . When R has a value of 

o 

0.025 cm the material tensile strength must be approximately 20 kilobars, the 
_--_es 1; imate d_dyn ami c^t ens i le_strength-of— the— target— material-,— 2 02 4-13-al : um-inum : .— 

As seen from Figure 20, there would be no failure in target No. 4; there would 
be a slight spall in target No. 3; a larger spall would occur in target No. 2; 
and the fracture strength would be exceeded throughout target No. 1 except near 
the rear surface. The momentum trapped between the spall fracture and the rear 
surface where the normal stress is always zero would undoubtedly cause complete 
penetration as observed in the experimental target. This model seems to give 
results that are in complete agreement with the experimental observations and 
the estimated dynamic tensile strength of the material. A summary of the ana- 
lytical and experimental results for these four aluminum targets is given in 
Figure 21.- The radial and tangential components of the stress waves in target 
No. 2 for various times are shown in Figure 22. 

The values of tensile stresses shown in the previous figures were for 
points along the target axis. Numerical values of the maximum stresses com- 
puted for points throughout these four targets are given in Tables II, III, 

IV, and V. The three causes of the tensile stresses are indicated on these 
tables. The heavy lines enclose all values exceeding the dynamic strength of 
20 kilobars. Data from these tables have been plotted in Figure 23. The lat- 
eral extent of the area in which the stresses exceed the critical tensile 
strength of the material can be determined by following the 20-kilobar contour. 

This analysis should predict failure of the material but does not describe 
it. If these targets had been some plastics the fractured area would probably 
be the entire area in which the fracture strength was exceeded (Reference 18). 
In aluminum targets, however, the fracture would appear as an internal crack 
following the portion of the 20-kilobar contour near the rear surface of the 
target (inner surface of the spacecraft mainwall). A part of the stress wave 
would be trapped between this fracture and rear surface producing a bulge or 
a rupture of the material. Also, this fracture would act as a new free surface 
for the tail of the pressure pulse, which, upon reflection from the fracture, 
might create additional internal cracks in the target (References 9 and 19). 


10. TENSILE STRESSES IN LAMINATED 
POLYCARBONATE-ALUMINUM TARGETS 


The relations given in Section 5 for the reflection and transmission of 
stress waves at an interface between two media were based upon a perfect bond 
between the two materials and no loss of energy at the interface. If there is 
perfect contact between the two materials but no bond, then compression waves 
will be transmitted but no tensile components can pass from the first material 
into the second. If there is not perfect contact between the materials or if 
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the dimensions of any voids between the materials are unknown, then a corre- 
lation between the computed stress and the observed fracture is impossible. 
Unfortunately, the two experimental targets Nos. 5 and 6 fit the last descrip- 
tion. This obviously results in a more efficient spacecraft mainwall than if 
they had been perfectly bonded unless such a bond is needed for structural 
reasons. It would have perhaps been even more efficient if the two materials 
had been separated by some porous filler material. Summers and Nysmith (Ref- 
erence 2) report that tests have indicated that low- density fillers such as 
polyurethane foam, glass wool, or aluminum honeycomb, can improve the effec- 
tiveness of a meteoroid bumper in protecting a space structure. DiBattista 
and Humes (Reference 20) have also shown that laminated mainwalls are superior 
to solid mainwalls in retarding spallation and penetration failure. 

A search of the literature failed to reveal any information about the 
momentum multiplication associated with high-velocity impact of polycarbonate 
targets. If it is assumed that the momentum multiplication factor for the 
impact of the polycarbonate first layer and the cavity radius are the same as 
for the debris-aluminum impacts , the computed stresses developed in the alu- 
minum would have been sufficient to rupture both the experimental targets. 
There was, however, no failure in target No. 6, but there was penetration 
failure of target No. 5. 

One must conclude that the maximum tensile stress in the second layer of 
target No. 6 did not exceed the dynamic strength of 20 kilobars. If the 
transmitted pressure wave amplitudes in the experimental targets were only 
20 percent of that transmitted through a perfectly bonded laminate, the com- 
puted stresses along the axis of the aluminum second layer would be approxi- 
mately as shown in Figures 24 and 25. The stress distributions throughout 
these targets are shown in Figure 26. 


11. EFFECTS OF BUMPER- MAINWALL SPACING 


The tensile stresses were computed for polycarbonate-aluminum targets 
where T^ = = 0.127 cm, for bumper spacings of 0.4, 0.8, and 1.2 cm. The 

results shown in Figure 27 indicate that for a spacing of 0.4 cm the target 
would be damaged extensively, perhaps perforated; that for a spacing of 0.8 cm 
there would probably be a slight internal fracture; and that there would be no 
damage for a spacing of 1.2 cm. As would be expected from the relations given 
in Equation (8), the maximum tensile stress varies approximately inversely as 
the square of the distance between the bumper and the mainwall. 


12. EFFECT OF VARYING LAMINATION THICKNESSES 
KEEPING THE MAINWALL WEIGHT CONSTANT 


Computations were made of the maximum tensile stresses in the aluminum 
laminate of a polycarbonate-aluminum target, keeping the total weight equal to 
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that of the 0.318-cm solid aluminum target of shot No. 3. The results are 
given in Figure 28 for three combinations of T and . It appears that the 

larger the ratio T^/^, the more resistant the mainwall would be to meteoroid 

impact. There is, of course, some limiting value to this ratio because of 
structural .requirements .. — - - - - - 


13. CONCLUSIONS 


All computed values and conclusions have been based on a 2017- T4 aluminum 
projectile and a 2024-T3 aluminum bumper. Applications to other materials 
should be made with caution. 

A bumper-protected polycarbonate-aluminum mainwall with no bonding between 
the two materials will increase the meteoroid- impact resistance of a structure. 
It appears that such a laminated mainwall will afford the same resistance to 
high-velocity impact as a solid aluminum structure weighing more than twice as 
much. Bonding of the two materials would result in greater damage and the use 
of a porous filler between the materials would probably result in even less 
structural damage.. A much greater resistance to penetration would be achieved 
if it were possible to make the outer laminate of aluminum and the inner one of 
polycarbonate . 

The conclusion that the amplitudes of the stress waves transmitted from 
the first to the second laminate of the experimental shots were approximately 
20 percent of that transmitted through a perfectly bonded laminate was based 
upon the very unlikely assumptions that the momentum multiplication factor for 
aluminum-polycarbonate impact is the same as for aluminum- aluminum impact and 
that the cavity radius, R Q , is the same for the two materials. 

It is suggested that a series of carefully controlled experiments be con- 
ducted for impacts in both bonded and unbonded laminates. The momentum multi- 
plication factor for debris-polycarbonate impact should be determined. 
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TABLE II. 


COMPUTED MAXIMUM TENSILE STRESS (KILOBARS) - SHOT NO. 1 


Y — 


(inches /cm) 


.X_(:inchesr/:cm) _ — 


0 . 0 / 0.0 


0.05/0.127 


0.10/0.254 


0.15/0.380 


0.015/0.038 

121.0* 

98. 8* 

49 . 0* 

0.020/0.051 

67.4* 

55.9* 

CD 

O 

CO 

0.025/0.0637 

44.1* 

37.5* 

22.4* 

0.030/0.0765 

32.2* 

27.9* 

17.7* 

0.035/0.089 


22.1* 

17.2*** 

0.040/0.101 

22 . 3*** 

23.1*** 

23.4*** 

0.045/0.114 

34.0*** 

32 . 7*** 

26.2*** 

0.050/0.127 

40 . 0*** 

35.4*** 

24.4*** 

0 .055/0.139 

33 . 0*** 

28.1*** 

17.2*** 

0.060/0.152 

14. 3*** 

11. 7*** 

6 . 74*** 


13 . 3*** 
10.4*** 
6.82** 
3.94** 


NOTES FOR TABLES II THROUGH V 


•'Tangential tensile stress near the target's front surface. 
**Radial tensile stress produced by the forward-moving wave. 
***Radial tensile stress caused by the reflected wave. 
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TABLE III. 

COMPUTED MAXIMUM TENSILE STRESS (KILOBARS) - SHOT NO. 2 


Y 

(inches/ 

cm) 


0.015/0.038 

0.020/0.0508 

0.025/0.0635 

0.030/0.076 

0.035/0.089 

0.040/0.104 

0.045/0.114 

0.050/0.127 

0.055/0.140 

0.060/0.153 

0.065/0.165 

0.070/0.178 

0.075/0.191 

0.080/0.204 



X 

( inches/cm) 


0.0/0 .0 

0.05/0.127 

0.10/0.254 

0.15/0.380 


121.3* 

67.3* 

41.5* 

32.2* 

20 . 6 * 
17.3* 
16.2** 
15.5** 
21.5*** 
29 . 6*** 
31.8*** 
24.4*** 
7.9*** 


98.5* 

56.0* 

37.4* 

27.9* 

22 . 0 * 


18.2* 


15.7** 


15.2** 


14.6*“' 


20.9*** 
27.7*** 
29 . 0*** 
21.5*** 


6.6 


49.0* 

30.8* 

22.4* 


17 . 7* 

14.6* 

13.0** 

12 . 6 ** 

12 . 2 ** 

14.6*** 

19 . 6*** 


22.6 

20.9 





8 . 3 * 
14.3*** 
15 . 0 *** 
13 . 9*** 
11 . 2 *** 

7 . 0 *** 
3.06*** 












TABLE IV. 

COMPUTED MAXIMUM TENSILE STRESS (KILOBARS) - SHOT NO. 3 


(inches/cm) 


_ 070X57 OT 0 38 
0 . 020/0 .0508 
0.025/0.0635 
0.030/0.076 
0.035/0.089 
0.040/0.104 
0.045/0.114 
0.050/0.127 
0.055/0.140 
0.060/0.153 
0.065/0.165 
0.070/0.178 
0.075/0.197 
0.080/0.204 
0.085/0.219 
0.090/0.228 
0.095/0.241 
0 . 100 / 0 . 254 ’ 
0.105/0.267 
0.110/0.279 
0.115/0.292 
0.120/0.304 


1 X (inches /cm) 

0 . 0/0.0 

. 

0.05/0.127 

0.10/0.254 

0.15/0.380 


121.3 
67.3* 
41.5* 
32.2* 
25.0* 
20 . 6 * 


17.3* 
16.2** 
15 . 5** 
14.8** 
14.1** 
13 . 4** 
12 . 7 ** 
12 . 1 ** 


11.5** 


11 . 1 ** 


10 . 6 ** 


12 . 1**M 


16 . 9'“'““ 


21.4 

20.6 


12.8 


98.5 

56.0 
37.4 
27.9 

22.0 


18.2 

15.7 

15.2 

14.6 

13.9 

13.2 

12.6 
12.0 

11.4 

10.9 

10.5 
10.0 

11.6 

16.3 


20.2 


49 . 0*~ 
30 . 8* 
22.4* 


17.7* 
14.6* 
13.0** 
12.5** 
12 . 2 ** 
11.7** 
11 . 2 ** 
10 . 8 ** 
10.3** 



8.35* 
8.1 ** 
7.86** 
7.66** 
7.45** 
7.25** 
7.15** 
6.93** 
6.75 ** 
6 . 68 ** 

10.4 *** 

12 . 4 *** 
12 . 4 *** 
10.1 *** 

5.45*** 










TABLE V. 


COMPUTED MAXIMUM TENSILE STRESS (KILOBARS) - SHOT NO. 4 


Y 

(inches/cm) 

X (inches /cm) 

0 . 0/0.0 

0.05/0.127 

0.10/0.254 

0.15/0.380 






0.015/0.038 

121.3* 

98.5* 

49.0* 


0.020/0.0508 

67.3* 

56.0* 

30.8* 


0.025/0.0635 

41.5* 

37.4* 

22.4* 


0.030/0.076 

32.2* 

27.9* 

17.7* 


0.035/0.089 

25.0* 

22.0* 

14.6* 


0.040/0.104 

20.6* 

18.2* 

13.0** 


0.045/0.114 

17.3* 

15.7** 

12.5** 


0.050/0.127 

16 . 2** 

15.2** 

12.2** 

8.35* 

0.055/0.140 

15.5** 

14.6** 

11.7** 

8.1 ** 

0.060/0.153 

14.8** 

13.9** 

11.2** 

7.86** 

0.065/0.165 

14.1** 

13.2** 

10.8** 

7.66** 

0.070/0.178 

13.4** 

12.6** 

10.3** 

7.45** 

0.075/0.197 

12.7** 

12. 0** 

9.9=** 

7.25** 

0.080/0.204 

12.1** 

11.4** 

9.6** 

7.15** 

0.085/0.219 

11.5** 

10.9** 

9.3** 

6 . 93** 

0.090/0.228 

11.1** 

10.5** 

8.95** 

6.75** 

0.095/0.241 

10.6** 

10.0** 

8.6 ** 

6.68** 

0.100/0.254 

10.1** 

9.65** 

8.3 ** 

6 . 47** 

0.105/0.267 

9.75** 

9.3 ** 

8.0 ** 

6.33** 

0.110/0.279 

9.35** 

8.9 ** 

7.8 ** 

6.2 ** 

0.115/0.292 

9.6 ** 

g % 3 AAA 

8.53*** 

6 . 06** 

0.120/0.304 

13.0 ** 

12.8 *** 

12.1 *** 


0.125/0.318 

17.5 ** 

17.0 *** 

15.0 *** 

12.1 *** 

0.130/0.330 

19.3 ** 

18.1 *** 

15.1 *** 

10.9 *** 

0.135/0.343 

14.9 

13.7 *** 

10.9 *** 

7.4 *** 

0.140/0.356 

4.86** 

4.4 AAA 

3.26*** 

2.33*** 






Figure 2. Debris Momentum 



Figure lj. Lemniscates of Bernoulli 
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Figure 5. Various Modifications of the Le ranis cate 
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0.01 0.02 0.05 0.1 0.2 0.5 1.0 


Example: Required a function that will reach a maximum 3tress in 0.5 micro- 

second and then decay to a value of 1/e of the maximum in 2.0 
microseconds. To use the chart, F = 0.10 and t D /t R = 4. Locate 
the point shown at t R /F = 5 and t D /t R =4. On the bottom scale 
read Faj_ = 0.055, giving = 0.55, and on the scale at left read 
A = 8.5, giving a ^ - 4.67. 

Figure 7. Chart for Determining Decay Constants 
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ANGLES OF TRANSMISSION a REFLECTION 



ANGLES OF TRANSMISSION 8 REFLECTION STRESS AMPLITUDE 
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Figure 16. Stress Waves at Interface and Free Surface 
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Figure 17. locations of Imaginary Cavities 
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radial tensile stress produced 



Figure 18. Possible Causes of Tensile Stresses 
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Figure 23. Continued 



Figure 23. Concluded 
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Figure 2h. Maximum Stress Along the Y-axis of Laminated Target - 

Shot No. 5 
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Figure 26. Maximum Tensile Stress in Laminated Targets 




Figure 26. Concluded 
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